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Symmetry-Protected Topological Phases at Finite Temperature 


O. Viyuela, A. Rivas and M.A. Martin-Delgado 

Departamento de Fisica Teorica I, Universidad Complutense, 28040 Madrid, Spain 

We have applied the recently developed theory of topological Uhlmann numbers to a representative 
model of a topological insulator in two dimensions, the Qi-Wu-Zhang model. We have found a stable 
symmetry-protected topological (SPT) phase under external thermal fluctuations in two-dimensions. 
A complete phase diagram for this model is computed as a function of temperature and coupling 
constants in the original Hamiltonian. It shows the appearance of large stable phases of matter 
with topological properties compatible with thermal fluctuations or external noise and the existence 
of critical lines separating abruptly trivial phases from topological phases. These novel critical 
temperatures represent thermal topological phase transitions. The initial part of the paper comprises 
a self-contained explanation of the Uhlmann geometric phase needed to understand the topological 
properties that it may acquire when applied to topological insulators and superconductors. 
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I. INTRODUCTION 

Topological insulators and superconductors w are novel phases of matter that cannot be described using the 
standard Gizburg-Landau theory for phase transitions. These novel phases are of great interest for applications 
in emerging fields such as spintronics, photonics or quantum computing. They have opened an area of new rich 
physical phenomenology, giving birth to materials that host robust metallic edge states and the possibility for realising 
Majorana fermions [Mill, as well as a zoo of new exotic particles known as ‘anyons ’ [HHH]. Actually, graphene was 
the first material proposed to realise a topological insulating phase using a mechanism of spin-orbital coupling between 
two graphene sheets [2]. Unfortunately, the gap due to this effect in carbon atoms is too small and the quantum spin 
Hall effect was finally achieved in HgTe quantum wells instead [IHIIIS]. Very recently, bosonic topological orders have 
been also proposed [2QH24] . 

Whereas these systems have been extensively studied for the idealized case of zero temperature, the common belief 
was that these exotic phases are gradually spoiled as the systems are heated up (25] [26] . In [25] , it was shown that 
certain topological properties like the edge states can be destroyed by a thermal bath even when the system-bath 
interaction respect the protecting symmetry of the topological insulator. It was not until the work in m that 
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a topological indicator with a dependence on the temperature of the system was introduced, that was related to 
the topological part of the quantum Hall conductivity and the degree of mixedness of the edge states due to the 
dissipative effect. Other attempts to understand the interplay between topological phases and finite temperature 
[28l|29], nonequilibrium dynamics [30] and dissipation-engeneering |3T]-[34] have also been made. 

The extension of the notion of symmetry-protected topological (SPT) order for fermions to situations where the 
system is not in a pure state, but in a quantum mixed state described by a density matrix has been developed very 
recently. The first realization of topological fermonic phases in the presence of thermal effects or external noise was 
proposed in Ref. [35]. Here a topological invariant remains quantized for a finite range of temperature/noise for the 
first time in emblematic examples of topological insulators and superconductors in ID. Then, new advances have been 
achieved for fermionic phases of matter with SPT order in two-dimensional systems [31137]. The novelty of these 
works relies on the introduction of two new observables: the Uhlmann phase for wire-shaped systems (ID) [35|, and 
the Uhlmann number for planar ones (2D) [37] as well as a variety of extensions [36], which are able to account for 
topological order even at non-zero temperature. Remarkably, these recent works report purely “thermal” topological 
phase transitions between distinct non-trivial topological regions solely driven by temperature, without losing their 
quantum character. It is also shown that there exists always a critical temperature above which the topological phase 
abruptly disappears. The possibility of having several topological thermal critical points has also been discovered as 
a novel effect m 

In this way, exotic topological quantum properties can survive even in the presence of dissipative effects and finite 
temperature, which are the natural conditions in real life experiments. This situation for fermion systems with SPT 
order is in sharp contrast to some results establishing the absence of stable topological phases subject to thermal 
effects, like for topological orders with spins [38ti42] . for SPT orders with spins m, or SPT orders with fermions 
under certain conditions [25]. Therefore, these new results open the way towards the characterization of fermion 
quantum phases of matter with topological phases protected by symmetry in thermal states, or more general density 
matrices. Note that other studies on thermal effects and related issues in topological systems have been recently 
carried out [43H52]. 

The aim of this paper is to give a complete and comprehensive derivation of the Uhlmann geometric phase as a 
tool to characterise symmetry protected topological orders at finite temperature, performing a detailed analysis of its 
intrinsic properties. Firstly, we will formally derive the expression for the Uhlmann phase as an extension of the Berry 
phase to the case of dealing with quantum mixed states described by density matrices. Later, we apply the formalism 
to new 2D models for topological insulators with time-reversal broken symmetry (TRB) that can be characterized 
using these Uhlmann geometric concepts by the so-called topological Uhlmann numbers m- 

The paper is organized as follows. Sec. [TT| introduces the concept of Uhlmann holonomies, giving a detailed derivation 
of the Uhlmann parallel transport condition and the geometric Uhlmann phase, by means of amplitudes and purifica¬ 
tions of a density matrix. In Sec. in the topological Uhlmann phase in ID is presented and applied to a generic ID 
two-band model for a topological insulator or superconductor. In Sec. |IV[ we derive the topological Uhlmann number 
for 2D topological systems, discussing some technical issues and using it to characterise an emblematic 2D topological 
insulator, namely the Qi-Wu-Zhang model [53] [54]. Sec.|v|is devoted to conclusions. 


II. UHLMANN HOLONOMIES 

Geometric phases have played an essential role in many quantum phenomena since its modern discovery by Berry 
[55] (see also Refs. |56l|57|). An emblematic example is the description of the quantum Hall effect using the theory 
of fiber bundles. More concretely, one can relate the transversal conductivity axy of the system with the first Chern 
number Ch as ISHIISS], 


a^y = Ch^, (1) 

where e is the electron charge and h is the Planck constant. 

The Chern numbers are topological invariants associated to characteristic classes of complex vector bundles when 
working on an oriented manifold of even dimension. The first Chern number Ch, which is defined for the 2D case, is 
just the integral of the Berry curvature over the two-dimensional parameter space, i.e. the Brillouin zone (BZ) in the 
context of condensed matter systems. In addition, this topological observable can also be related to the presence of 
conducting edge states when the system has open boundary conditions. 

The standard way to observe the properties and implications that the Berry phase has in solids (Zak phase [60]), 
has been through indirect measurements like the transverse conductivity, or side effects like the presence of conducting 
edge states at the boundary. However, very recently the first direct experimental measurement of a Berry phase in 
a cold atoms system has been achieved [6TH63] . In particular, m reports the measurement of the Zak phase in a 
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one-dimensional optical lattice simulating the different phases of polyacetylene, characterizing the topology of the 
system and observing the quantization of this phase. This paves the way towards the extension of geometric phases 
to study topological properties beyond the realm of condensed-matter systems [64f[67] . Moreover, there are recent 
experiments on this context involving architectures of superconducting qubits [68] [69l . 

Although these studies provide a successful picture for quantum systems in pure states, typically the ground state, 
very little is known about the fate of those topological phases of matter when the system is in a mixed quantum state 
represented by a density matrix. 

A fundamental problem in the theory and applications of geometrical phases is the extension from pure quantum 
states (Berry) to mixed quantum states described by density matrices. Uhlmann was first to mathematically address 
this issue m and to provide a satisfactory solution mHni. For more than a decade, there has been a renewed interest 
in studying geometric phases for mixed states and under dissipative evolutions from the point of view of quantum 
information ffStiSO] . In the rest of the section, we will give a detailed derivation of the Uhlmann geometric phase 
and the definition of Uhlmann parallel transport. These two key ideas will allow us to construct later two topological 
observables to distinguish different topological phases in condensed matter systems both in ID and 2D. 


A. Amplitudes and Purifications 

One of the most fundamental problems when applying geometric concepts to the set Q of density matrices is the 
nonlinear character of it. That is, consider some scalars Ai and A 2 , for all pairs of density matrices pi and p 2 , the 
linear combination Aipi + X 2 P 2 is not another density matrix unless Ai, A 2 > 0 and Ai + A 2 = 1. In the case of pure 
states p = |'0){'0|, this is not a major issue as we can immediately introduce a Hilbert space formed by the vectors \ 2 p) 
so that the quantum state p = |'0)('0| is the projector on the subspace spanned by jV^). Then the geometric properties 
of the Hilbert space can be applied to the vectors jV^), but having in mind that \%Ij) and describe the same 

physical state p = |'0){'0|. Namely, to work with a Hilbert space we have to pay the price of introducing an additional 
degree of freedom, the phase of the vectors, such that quantum pure states as seen in the Hilbert space are defined 
modulo a U(l) phase factor. 

The key idea is to realize that an analogous construction was possible for any general density matrix. Specifically, 
consider a density matrix p and define its amplitude as any of the matrices w such that 

p = ww^. (2) 

The amplitudes form a Hilbert space 1-Lw with the Hilbert-Schmidt product {wi^W 2 ) = Tr(re|re 2 ). From Eq. Q, we 
see that there is a U(n)-gauge freedom in the choice of the amplitude (here n denotes the dimension of the space), 
since w and wU^ where U is a unitary operator, are amplitudes of the same state. 

Note the parallelism between the usual U(l)-gauge freedom for pure states, and the U(n)-gauge freedom for the 
amplitudes, see Fig. In the case of pure states, one defines a projection map it : H ^ Q: 

7r(|^)) = |^)(V’|, (3) 

so that lip) and e'‘^\tp) represent the same physical state, i.e. the same density matrix given by the projection map tt. 
In the case of amplitudes, the projection map is analogously defined, tt : Hw ^ Q' 

7r(w) = wwK (4) 

Thus, the usual gauge freedom for pure states can be seen as a particular case of the amplitude U(n)-gauge freedom 
for general states. 

So far the concept of amplitude seems just a mathematical abstraction, however it also enjoys a nice physical 
meaning. To see it, consider what is called the “square root” section of the Uhlmann fiber bundle, which is nothing 
but the parametrization w = for amplitudes of some density matrix p (note that this is also called polar 

decomposition of re, see below). Then by using the spectral decomposition, p = we have that 
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Now we define the following isomorphism between the spaces 1-Lw and H (8) H, 
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U(l) U(n) 




FIG. 1: Comparison between the usual fiber bundle of pure states (a), and the fiber bundle proposed by Uhlmann (b) for any 
density matrix through the concept of amplitude. 


where denotes the transposition of U taken with respect to the eigenbasis of p. Then, the property p = is 
now written as 


p = Tr2(|«;)(u;|). (7) 

Here, Tr 2 denotes the partial trace over the second Hilbert space oil-L^l-L. In other words, any amplitude w of some 
density matrix p can be seen as a pure state \w) of the enlarged space H ^ 1-L with partial trace equals to p, Eq. 
0. Thus, it is said that \w) is a purification of p. Thanks to this isomorphism between amplitudes and purifications 
the whole geometric approach by Uhlmann can be written in terms of purifications which are objects with a clear 
physical meaning. However, we shall keep on developing the theory and explanations using amplitudes and refer to 
the purification interpretation when needed. 


B. Parallel Amplitudes 


The concept of parallel amplitude is essential to define the Uhlmann parallel transport condition and the holonomy. 
Given a pair of states pi and p 2 , we define two of their amplitudes wi and W 2 to be parallel if they minimize the 
Hilbert space distance in i-e. wi \\ W 2 if 


ll'^i - '^ 2 |P = min \\wi - re 2 |p, with pi = wiwl, p 2 = W 2 wl. 

W\,W2 


We may work further on this condition and establish some properties of such wi and W 2 . 

min \\wi — W 2 \\‘^ = min {wi — W 2 ^Wi — W 2 ) 

Wi,W2 Wi,W2 

= min Tr ( w\wi + W 2 W 2 — W 1 W 2 — w^wi ) 
Wi,W2 V / 

= Tr(pi) + Tr(p 2 ) — max Tt{w\w 2 + w^wi) 

Wi,W2 

= 2 — 2 max Re Tt{wIw 2 ) . 

Wi,W2 L 


(8) 


(9) 


Since Re(x) < |x|, it is clear that the condition for the maximum is to chose rci, W 2 ^ such that such w\w 2 is self-adjoint 
and positive definite 


w\w2 = W 2 W 1 > 0 . ( 10 ) 

More explicit expressions may be obtained with the polar decomposition theorem. This theorem states that any 
operator A (consider it to be finite dimensional for simplicity) can be decomposed as A = where |^| = VAA^ 
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and is a unitary operator [SUIH]. Furthermore for any unitary operator U the following inequality holds 
Re[Tr(A[/)] < \Tv{AU)\ = \Tr{\A\UUA)\ 

= Tr I < ^iTr\A\) [if (u^^W\A\UUa)] = Tr\A\ (11) 

where in the penultimate step we have used the Cauchy-Schwarz inequality \Tt{A^B)\‘^ < Tt{A^A)Tt{B^B). Since 
the equality is reached for U = we can write 

maxRe[Tr(A[/)] = Tt\A\. (12) 


Now, the polar decomposition theorem on the amplitudes wi = ^/^Ul and W 2 = y/^U 2 jointly with Eq. (12) lead to 

Tt{wIw 2 ) = max Re TT{Ul^/^^/^U 2 ) 

J Ui,U2 L 

= max Re [Tr{^^U)] = Tr\y^^\ = Tr^J vWp2vW- (13) 


max Re 

Wl ,W2 


Therefore, the Hilbert space distance between two parallel amplitudes wi and W 2 is equal to the Bures distance 
between pi and p 2 , 


Iki - 1^2 IP = 2 - 2 Tr^for Wl II ^2. 


(14) 


Finally, note that the equality in (13) is obtained for U = this is the adjoint of the unitary operator of the 

polar decomposition of ^ unitary operator is not uniquely defined unless y/piy/^ 

is full rank (and so y^, and \^/^^/^\). In such a case, is given by the equation 


(15) 


Since from Eq. (13) we have U = U 2 UI = by using the unitarity of U and Eq. (15) we arrive 


at 


(16) 


Hence, Eq. (16) gives us the relation that the unitaries Ui and U 2 must satisfy in order to minimize the distance 
between the associated amplitudes wi and 1 ^ 2 , i.e. in order to have wi || W 2 . 


C. Uhlmann Parallel Transport and Connection Form 

Before introducing the Uhlmann holonomy, we need to explain some of its ingredients: the Uhlmann parallel 
transport condition and the connection form. Consider a trajectory in the density matrix space Q, pr- This induces a 
trajectory in the space via the amplitudes Wr^ pr = Wrwl. Since for a point pr of Q there are several amplitudes 
Wr, a unique trajectory on Q induces multiple trajectories in Bw depending on which Wr are chosen. To select one of 
them, we define a parallel transport condition, i.e. for a point pr with some amplitude Wr^ we consider the amplitude 
Wr-\-dr of the next point of the trajectory, Pr+dr, to be the closest to Wr according to the Hilbert space distance, so 
Wr II Wr-\-dr lu the scusc of the previous section, see Fig. Such a path in Bw is known as a parallel lift of path pr 
in Q. This is analogous to the Berry parallel transportation, where the U(l) phases of the consecutive state vectors 
\f^r) and \ f^r+dr) are chosen such that the distance between them ^/(('^rl — (V^r+drDdV^r) — IV^r+dr)) is minimal. 

Let prQ be the initial state of the trajectory, under the square root section (or equivalently by the polar decomposition 
theorem) we have Wr^ = s/pr^Uro- For the sake of clarity we shall work for the moment on a gauge such that Ur^ = L, 
and we shall refer to this gauge as the “identity” gauge. At some point r of the trajectory we have Wr = y^Ur, 
where Vr is a unitary (Uhlmann holonomy) constructed from the initial condition WrQ = sjpro by applying the 
parallel transport rule along the trajectory. The infinitesimal generator A{r) of this unitary transformation Vr is 
called connection form of the trajectory. Thus, if r{t)\l^Q is a parametrization of the trajectory [r(0) = ro], Vr is 
fully determined by the differential problem 

\V(0) = 


( 17 ) 
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FIG. 2: Uhlmann parallel transport condition. The amplitude Wr+dr is the closest amplitude to Wr among the ones in the 
hber at r + dr. 


whose formal solution can be written independently of the parametrization as 


X. ^ r A(r')dr' 


(18) 


where V is the path ordering operator. 

As we shall see, the Uhlmann connection A is not uniquely defined unless the density matrix pr is full rank. In 
order to compute A in such a case, let us rewrite Eq. in terms of differentials 


di/yt = 


(19) 


where A = ^ • Aidvi is understood as a differential one-form. Note that in Eq. 0, A[r{t)] = Then, 

consider some full rank density matrix pi = p, and let displace it by an infinitesimal translation p 2 = P d- dp. 
The basic condition to compute A is that the corresponding infinitesimally closed amplitudes wi = s/pV and W 2 = 
\/pd- dp{V + dU) fulfill the parallel transport condition (p^. Thus, Eq. (16) leads to 


{V + (W)v‘'= {-/p + A-Zp) '^y/p-^J^{p + Ap)^ 


( 20 ) 


where we have used that for an infinitesimal translation dp, p ^ dp = ^p+d^p. Eollowing Hiibner [73], we introduce 
an auxiliary real parameter s multiplying every differential: 


{V + sdV)V = (y^ + sdy^) ^ Vp{p + sdp)y/p. 


( 21 ) 


The aim of this trick is to expand the above equation around 5 = 0 and keep only the first order in 5. In order words, 
we have that 


1 + sdi/yt = 1 + 


d 

ds 


(VP + sdv^) ^^/p~^^./p{p +sdp)y/p 


s=0 


therefore 


A = dVV^ = 


-(v^ + sdv^)-i 


Vp + p 


-1 


s=0 


■^\JVP{P+Sdp)y^ 


( 22 ) 


(23) 


s=0 


To give a closed expression for A, we determine its matrix elements with respect to the eigenbasis {{'ipj)} of p, which 
we assume to be diagonalized with eigenvalues {pj}. Eor the first term in (23) we have 


1 -1 


{^ + sd^/p) ^(1 + svT^d^) ^ 

= (1 + s^/p^d^/p)~'^^/p-^y/p = 1 - s^/p-^d^/p, 


( 24 ) 
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where in the last step the inverse has been taken at first order in s. Hence 


-(^p + sd^p)-^ 


y/p= 


s=0 


and in the eigenbasis of p, 


- = — l={^i\d^\^pj). 

yPi 


(25) 


(26) 


The second term of (23) is a bit more involved. Define K{s) := we need to compute p 

By differentiating the square, 


{[K{s)fy = K'{s)K{s)+K{s)K'{s) = .fpdp^p, 

as i^(0) = p we have 

K'{d)p + pK'{d) = ^pdp^p. 

After taking matrix elements in the eigenbasis of p, 

{Pi + Pj){'>Pi\K' {d)\'ipj) = ^/PiPj{yi\dp\yj). 

Thus, the matrix elements of p~^K'{d) are 

mp-^K'myy = i(t/,,|i^'(o)i^,-) = 

Now, since dp = <^{-s/p-s/p) = (dy/p)^/^ + Y/p(dy^), we obtain 


-Ii^/^OM / \ _ V^(\/^+vW) 


y/Pi{Pi+Pj) 


{yi\d^\yj) 


and finally the matrix elements of the connection form. 


(^il^lV’i) = { —]= + {'tpi\d^\'tpj) = ^^(^iid^iV’i). 


y/Fi{Pi+Pj) 

This can be rewritten in terms of a commutator. 


{pi+pj) 






ipi +Pi) 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


Note that this connection form has only zeroes on its diagonal and is skew-adjoint, so that the Uhlmann connection 
is special unitary. 

For the particular case of 2 x 2 density matrices, because of the normalization condition Pi +P 2 = 1, Eq. (33) can 
be written as A = [dy/p, ^p]. This formula was also obtained by Dittmann and Rudolph [84]. Moreover, expanding 
in terms of Pauli matrices y/p = uq + a • cr, we may write A = 2i(da x a) • cr [73] . 

The result ( [3^ has been obtained under the requirement that p is full rank. If such is not the case, for wi = Wr the 
formula (14) can be fulfilled for several W 2 = Wr-ydr and there is not a unique solution and so a unique connection. 
One of these solutions can be obtained from continuity if we have a family of full rank states px which approaches 
the desired non-full rank state as A ^ 0. Then a connection can be obtained with the formula (33) for p^, taking the 
limit A ^ 0 at the end of the computation. We will come back to this point at the end of the following section. 

Once a connection form A is defined, we can construct its associate curvature form via the standard relation. 


F = dA + A A A, 


(34) 


where A denotes the exterior product and the action of d is understood as the exterior derivative [85] . For the sake of 
clarity, let ns write this equation in components. Expanding A = A^dr^, we find the usual expression of the field 
strength. 


Fij — diAj — djAi + [A^, Aj], 

ct the Chern classes of the 
way to obtain topological invariants. We shall come back to this point in Sec |IV| 4. 


(35) 


where F = ^ -Pijdr, A dr^. 


2 ^i,j ^ 

This curvature form F allows us to construct the Chern classes of the Uhlmann fiber bundle, which is the standard 
















D. Gauge Transformations 


In the previous subsection we obtained an expression for the Uhlmann connection in the “identity” gauge, i.e. for 
the trajectory pr, we take initially = ^/Pr{o) construct the next point of the trajectory in Hw applying the 

Uhlmann parallel transport rule, such that Wr = ^/^Vr where Vr is the Uhlmann holonomy of the trajectory. Now 
consider a change of gauge by left-multiplying the amplitudes by some r-dependent unitary Ur^ Wr^ = ^Jpro Uro , so 
that Wr = y/^UrVu^n whcrc Vu^r is the Uhlmann holonomy in the new gauge. To compute Vu^r^ we first obtain the 
connection form in the new gauge. 

Similarly to the previous section we have to impose that two infinitesimally closed amplitudes are parallel. Under 
an infinitesimal displacement wi = y/pUVu changes to W 2 = \/p + dp{U + dU){Vu + dVjy), so that according to Eq. 
wi\\w2 if 


{U + dU){Vu + dVu)V^U^ = vU'\/x/p(p + dp)^. 


(36) 


From Eq. (20) and subsequent ones, the right hand side is simply 1 + where At denotes the connection in the 
“identity” gauge. Moreover, the left hand side at first order in differentials is 1 + dUU^ + UdVuV^U^, so we obtain. 


dUU^ + UdVuVj\u^ = Ai. 


(37) 


Since by definition, Au = dVjjVX is the connection form in the new gauge, we find the transformation rule 


Au = f/Uif/ - f/Mf/. (38) 

In general, under a change of gauge U 2 = UiU {w ^ wU)^ we straightforwardly obtain the transformation law 

Au2 = U'’Au,U -U^dU. 

Now we can calculate the holonomy in the new gauge, 

X. ^ r Au(r')dr' 

To that aim, consider the parametrization r{t)\l^Q and compute the t-derivative of U[r{t)]Vu^r wU^rm- 
I {U[rmVu,r(t)U^[rm} = + U[r{t)]^^^U[r{0)] 

= + U[r{t)]Au[r{t)]Vu,rit)U[rm 

= At[r{t)]U[r{t)]Vu,rit)U[rm, 


(39) 

(40) 

(41) 


(42) 

(43) 


where we have used the definition Eq. (40) and introduced the gauge transformation Eq. (38). Since 
U[r{t)]Vu^r{t)U^r{0)]\t=o = 1; {U[r{t)]Vu^rit)^[i^{^)]} satisfies the same differential problem as Vt^rit)^ Eq. ( pj| ), 
and the uniqueness of the solution implies that: 


= U[rit)]Vu,r(t)U^[rm ^ Vu,r = UHr)Vt,rUiro). 


(44) 


Therefore, the Uhlmann holonomy of the trajectory is clearly gauge covariant as expected. This is a fundamental 
requirement in order to define gauge invariant quantities out of this construction. 


E. Uhlmann Geometric phase 

Once the concept of Uhlmann holonomy has been presented, it is interesting and convenient to define gauge 
invariant quantities associated to it. These are constructed out of the relative phases of the amplitudes that have 
been parallelly transported. Actually, this is the most important part of the analysis, since we want to relate the 
mathematical concepts to physical observables that are represented by gauge invariants in the theory. 

First of all, consider some trajectory on Q, with initial and final points po and pi respectively. For some initial 
amplitude icq, po = wqWq^ a parallel lift leads to some amplitude rci, pi = wiw\. The Uhlmann geometric phase is 
defined as 


$u := arg(«;o,wi) = argTr(«;^u;i) 


(45) 
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U(l) U(n) 




FIG. 3: Comparison between Berry and Uhlmann Holonomies. In the approach of Berry, valid for pure states (a), after a closed 
loop on the space Q, the hnal state vector differs by a phase with respect to the initial state vector. In the case of Uhlmann, 
formulated for density matrices (b), the difference between initial and hnal amplitude is, in general, the unitary matrix given 
by ^. 


which is a gauge-independent quantity. Indeed, in the identity gauge, wq = and wi = we have 

$u = argTr(y^y^r). 


(46) 


(44), V U;VUo, so $u ^ 


Under a gauge transformation wq wqUq, wi wiUi and according to Eq. 
argTr(I7dy^y^[/i[/|UUo) = ^u- Therefore ff>u is gauge invariant. 

To better visualize the geometric meaning of ff>u, let us consider the trajectory on Q to be closed. Then po = pi = p 
and we have 


ff>u = argTr(pU) = argTr (^pVe^ . (47) 

In the Uhlmann approach, after a closed trajectory in Q the amplitude gains an extra geometric U(n)—phase, 
Ve^^. This is the equivalent to the Berry geometric U(l)—phase in this non-abelian approach, see Fig. The 
equation defines a real gauge invariant quantity extracting the non-abelian U(n)—phase which can be understood 
as the relative phase between wq and wi. Basically, this is because by using the amplitude-purification isomorphism 
Eq. (|^, can be equivalently written as 

ff>u = arg(u;o|u;i). (48) 

Actually, this interpretation has suggested some schemes to measure Uhlmann geometric phases m 

Another important concept needs to be introduced in order to understand the properties of the Uhlmann phase. 
This is the difference between an absolute holonomy and a pointed holonomy. The mapping 

Mu : w —U, (49) 

is a so-called pointed holonomy. This means, that even if the trajectory in Q is closed, in general the holonomy 
depends on the initial point of the path. Nonetheless, we have identified cases in which the pointed holonomy reduces 
to an absolute holonomy becoming independent of the initial point |35l |37] . In addition, it is always true that the 
mapping 


Mu : w —^ spect(U), (50) 

where spect(U) stands for the spectrum of the unitary U, is an absolute holonomy. 

To conclude this section, consider the case of a pure state trajectory pr = |'0r-)('0r|- As we have already pointed, 
the Uhlmann connection in this case is not uniquely defined because the state pr is not full rank. There are two 
approaches to overcome this difficulty: 


• As aforementioned, a first approach to compute the Uhlmann phase factor in this situation is to consider the 
connection computed via (33) for a full rank matrix trajectory p^ such that Mmx^o Pr — l'0r)('0r|- Thus, 
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the Uhlmann phase factor is obtained by formula (47) with V = 
when 


This is a typical approach 


is seen as the ground state of some Hamiltonian and we want to study temperature effects. Then, 
the family of full rank states will be the Gibbs states which approach the ground state in the zero temperature 
limit. 


• An alternative approach is to look for another solution to the optimization problem (14). For example, consider 
two pure states, IV^i) and |'02): the optimization problem (14) can be easily solved in a gauge with amplitudes 
of the form of wi = |'0i)(a|e^'^b '^2 = where is some arbitrary r-independent vector, and the 

phases are chosen such that (^i — ^ 2 ) = arg('0i|'02). This can be directly checked in (14) and also derived from 


condition (10). Thus, in this gauge we can satisfy the parallel transport rule just by varying a U(l) factor, e^"^, 
extracted from the general U(n) gauge freedom. More concretely, we associate the amplitudes Wr = \^r) 
with the trajectory pr = |'0r)('0r|- Then we have 


vl,{wr + dwr) = |a)(a| + e |a)('0r-|d 

= [1 + {^r\di^r) + id0(r)] |a)(a|. 


(51) 


where due to normalization, {iprld'ipr) is purely imaginary. The parallel transport law in the form of (10), with 
wi = Wr and W 2 = Wr dWr leads to the condition 


d(/)(r) = i{'ipr\d'ipr) • 


(52) 


If we denote I'lpro){'^ro\ the initial and final state, we have wq = l^pro )and the parallel transported 
final amplitude will be wi = \'ipro )where is the accumulated phase after the parallel 

translation. Thus, by using Eq. ( [5^ the Uhlmann phase factor 4>u = argTr(reQrei) yields 


4>u = argTr e \a) {^proli^ro) = i (p {^prld^^r) niod 27r 


(53) 


which is nothing but the usual Berry geometric phase. 

Very importantly, since the Uhlmann phase is independent of the gauge, the two approaches must provide the same 
result. This is going to be crucial in the next section in order to obtain consistent results when taking the limits of 
zero temperature and pure states. 


III. DENSITY-MATRIX TOPOLOGICAL INSULATORS AND SUPERCONDUCTORS: TOPOLOGICAL 

UHLMANN PHASE 

In this section we will describe how the Uhlmann geometric phase acquires a topological character in certain ID 
condensed matter systems. The problem of combining the study of topological phases with the fact of having mixed 
states and finite temperature is a very relevant and had remained vaguely explored for a long time. In a set of recent 
papers I251E71I351I37], we tried to deal with this issue and proposed some methods to generalize the notion of SPT 
order to finite temperature and to systems with dissipation. 

The most characteristic feature of a geometric phase, as the Berry phase, is the fact that it does not depend on 
properties such as the gauge or the speed at which a given path is traversed but only on the geometry of the path. 
In the case of closed paths, this is manifested by the fact that the geometric phase only depends on the area enclosed 
by the trajectory on the manifold defined as parameter space. 

There are cases, however, where the path covered during the holonomy coincides with one of the topologically non¬ 
trivial paths of a given manifold. In particular, in condensed matter, we know that due to translational invariance, 
the crystalline momentum k defines the so-called Brillouin Zone (BZ) and characterizes the eigenstates of the system. 
Then the phase difference picked up by the system when parallelly transported from an initial momentum ko to a 
final one ko + G, where G is a reciprocal lattice vector, is gauge invariant and a geometric phase too. This was proven 
by Zak m and the arguments work both in ID and in higher dimensional systems. This particular Berry phase is 
also called Zak phase. Actually, in ID systems, Zak [60] showed that the Zak phase is quantized and equal to 0 or tt 
in the presence of inversion symmetry. 

On the other hand, the Uhlmann phase was initially proposed as a generalization of the Berry phase to density 
matrices, but it can also be applied to condensed matter systems in the spirit of the Zak phase. Actually, this is 
very fundamental since it deals with the problem of generalizing the notion of topology to systems with dissipation 
and finite temperature. Indeed, we studied this problem in detail [35] and found that when applied to representative 
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models for topological insulators and superconductors, the Uhlmann phase acquired a topological character and is 
also quantized to 0 or tt depending on whether the system is in a topological phase or not. 

Specifically, we consider ID Hamiltonians with two bands in momentum space. We denote by a and b two different 
types of fermionic band operators respectively, in such a way that in the spinor representation of ipk = for 

insulators, while for superconductors, 'ipk = (c/c, is the Nambu spinor for paired fermions with opposite crystalline 
momentum [9]. We take lattice spacing a = 1 throughout the text. The Hamiltonian of these systems can be written 
as 


H = J24HkA, (54) 

k 

where is a quadratic form on ^ and Hk is a 2 x 2 matrix. We aim at studying the properties of this Hamiltonian 
under a displacement in k along a circuit which encloses the whole Brillouin zone, which will be the circle S^. 

Note that Hk can be decomposed in the following form 


Hk = f{k)l + • cr. 


(55) 


Here, a = (cr^,, cr^, (7;^) are the Pauli matrices, Ak corresponds to the gap of Hk and f{k) denotes some function of k. 
The unit vector rik = (sin 6> cos 0, sin 6> sin 0, cos 6>) is called winding vector^ where 0 and 0 are /c—dependent spherical 
coordinates. 

For a thermalization process preserving particle number, with the Fermi energy set in the middle of the gap, the 
connection in Eq. (33) associated to this equilibrium (thermal) state, = e“^^^/Tr (e“^^^) where T = 1/[3 denotes 
temperature, is 


A = mi 2 (u- \dku\) \u^-) {u^ \ dk + h.c. 


(56) 


where := 1 — sech ('^) 


Moreover, making use of the discrete symmetries uniiHzi that restrict the movement of the winding vector rik as 
a function of k, only two components and n^, with i ^ j are non-zero. In order to obtain a simplified expression 
for the Uhlmann phase, we fix the gauge of the eigenvectors of Hk in Eq. (55) in such a way that the off-diagonal 
overlap {u^ \dkU^) and the winding vector components are related in this manner 


in’! \dkul) 


dknj 

2ni 


(57) 


Note that as Eq. 
Uhlmann phase. 

Thus, using Eq. 
expression 


(47) is gauge invariant, we can fix the most convenient gauge to perform the 


(56) with this particular gauge choice and simplifying Eq. (47), we obtain 


calculation of the 
a more simplified 


4»u = arg 


cos(7rcc;i) cos 



(58) 


where ooi := ^ f winding number of Uk [35], a topological indicator characterizing the mapping 

5^, from k —n/^. 

In the case of being in a trivial phase with uji = 0, the Uhlmann phase is always zero. However, for nontrivial 
topological regions uji = ±1, we can obtain either 4>u = 0 or 4>u = tt depending on the value of T. Thus, the 
topological order as accounted by 4>u can be also non-trivial for a range of temperatures T < Tc. In [35| we computed 
4>u finite temperature for the three representative models models of topological insulators and superconductors, 
showing this critical effect with temperature represented in different phase diagrams. 

On the other hand, the appearance of a critical temperature Tc in the Uhlmann phase can be interpreted geomet¬ 
rically in a pictorial and intuitive way. Eirst of all, for the sake of illustration, we can visualize the amplitudes (or 
purifications) as arrows with fixed length. The angle between two arrows will be given by the relative phase between 
the corresponding amplitudes. In Eig. we sketch different situations where the amplitudes (arrows) are transported 
according to the Uhlmann’s parallel condition along a closed loop that covers the whole Brilluoin zone . Here we 
present a summary of the different behaviors: 
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ko 

kf 



ko 

kf 


0 

27r 


B.Z 

(b) 


FIG. 4: Behavior of the Uhlmann phase during the parallel transport from /co = 0 to /c/ = 27r. We identify two situations: 
(a) Non-trivial topology highlighted by the presence of a topological kink at kc. (b) Trivial order where the topological kink 
disappears. A detailed explanation of the figure can be found in the main text. 


1. At T = 0 two things may happen. If the system is in a trivial phase, then the arrow is transported pointing 

always the same direction, which means that the relative phase between the amplitudes does not change. The 
initial and final arrows are parallel as in Fig. lilb), SO that = 0. In the case that the system is in 

a topological phase, the Uhlmann phase remains 0 (the arrows are parallel) up to a certain critical point /Cc, 
when the direction of the arrows is reversed and hence the relative phase between amplitudes is equal to tt. 
This constitutes a topological kink. From there, the arrow does not flip anymore until the final point of the path 
kf = 27r is reached, as in Fig. [^a). Thus, the phase = tt. 

2. At T 7 ^ 0, we can identify three different possibilities. If the system was in a trivial phase at T = 0 then 

= O 5 VT, as we have already proven, being again in the case of Fig. [^b). If at T = 0 the system was in 
a topological phase, then at finite temperature there are two options. Whenever T < Tc, the situation in Fig. 
Qa) is reproduced, the Uhlmann phase remains 0 up to the critical point which is a function of T. From that 
point, the arrow is reversed and the Uhlmann phase is equal to tt until the close path is completed. An increase 
of temperature shifts the position of kc towards the end of the Brillouin zone. However, if T > Tc, then we 
recover again Fig. [^b) with no phase change = 0- The position of kc has already reached the end of the 
Brillouin zone and the topological kink can no longer occur along the path between ko = 0 and kf = 27 r. 

To conclude, we have derived a connection between the existence of a critical temperature Tc and the presence of 
a topological kink structure [ 88 ] in the Uhlmann holonomy. Actually, the same relation can be stablished with the 
parameters of the Hamiltonian that drives the system out of the topological phase. This is because the Uhlmann 
phase places on equal footing T and these coupling constants. The reason to call this kink topological is actually 
well established. The critical momentum kc divides the holonomy in two different topological sectors according to 
the value of the final point of the path kf. The first sector covers the region with kf < kc and <Fu = 0, and the 
second sector starts when kf > kc and = tt. It is clear that these two different sections cannot be smoothly 
connected without crossing the singular point kf = kc- In the trivial topological regime = 0 however, we can 
associate every point in the sector with the same class. To be more specific, is invariant under changes of the final 
point of the holonomy kf e (0,/Cc). The same happens within the other non-trivial topological sector with <I>u = tt 
once we consider kf G (/Cc,27r). The effect of the temperature in the Uhlmann parallel transport can be understood 
as a displacement of the topological kink along /c-space. When the temperature reaches the critical value Tc, the 
critical momentum kc has been pushed out to the edge of the Brillouin zone and the topological kink disappears. 
Complementary, the critical temperature Tc can be viewed as the maximum amount of noise/disorder such that the 
Uhlmann holonomy along the Brilluoin zone displays a topological kink structure. 


IV. 2D UHLMANN TOPOLOGICAL INVARIANTS 

In the previous section, we have shown that in ID condensed matter systems, the Uhlmann geometric phase 
acquires a topological character when applied to topological insulators and superconductors. Actually, the relation 
of topological phases for pure states with invariants of SPT order in more than ID has been extensively studied, in 
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particular see nano]- Nonetheless, the concept of a geometric phase or a topological phase is not at all restricted to 
one-dimension and can be defined in higher-dimensional spaces. 

The problem of generalizing the Uhlmann phase formalism for 2D cases is not straightforward at all [36l [37] . This 
has to do with the triviality of the Chern classes associated with the Uhlmann connection of non-singular density 
matrices, as was first proven in Ref. [35]. Different but equivalent ways to solve this problem were proposed in Ref. 
[361137]. In Ref. [36] an analysis of the spectrum of the whole holonomy matrix was carried out, while in Ref. m 
a unique topological invariant called the Uhlmann number was first proposed to characterize the topological order of 
two-band topological insulators and superconductors at finite temperature in 2D. 

In this section we review the theory of topological phases for mixed states in 2D, and give a detailed explanation of 
a particularly useful way to define topological invariants inside the Uhlmann formalism. As a concrete example, we 
will characterize the topological phase diagram of the QHZ model for a topological insulator at finite temperature. 


A. Triviality of the Uhlmann Chern classes 

At zero temperature, there is a standard way to define topological invariants in two-dimensional symmetry protected 
topological order systems out of Chern numbers. In this case, a fiber bundle is formed out from the pure states 
eigenvectors \uk) corresponding to the lower energy band, and the fiber is linked to the U(I) gauge freedom that each 
eigenvector has. As aforementioned, the first Chern number associated to the first Chern class of fiber bundle can 
be related to the transverse or quantum Hall conductivity of the system by the so-called TKNN formula Eq. 0, 
[58] . Moreover, the presence of conducting edge states in topological insulators, or Majorana fermions in a topological 
superconductor, can be also related to the fact of having a topological non-trivial structure in the system. 

We shall consider a time-reversal broken two-band system with the Fermi energy between both bands. The Chern 
number is given by 


Jbz 


Ch := — / d?kF,y{k), F,y{k) := 


dyAl 


\k), 


(59) 


where BZ stands for Brillouin zone, A^{k) = i{uk\djUk) is the Berry connection and \uk) is the eigenvector corre¬ 
sponding to the lower energy band. This number is a topological invariant which only takes on integer values. While 
this kind of constructions can be extended to higher dimensional systems or systems with time-reversal symmetry 
[9], when attempting the generalization to density matrices via the Uhlmann connection, one finds the following 
fundamental obstruction. 

A natural way to generalize the Chern number to arbitrary density matrices would be to consider the first Chern 
class associated to the Uhlmann curvature, which is constructed from the Uhlmann connection A via the standard 


formula for the non-Abelian case, Eq.(35), 


^xy — 9xAy 


Then, using the theory of characteristics classes 
be given by 


dyAx-\-[Ax^Ay]. (60) 

EB, the (first) Chern number of the Uhlmann curvature would 




(61) 


There is however, a problem with this construction, as Chu turns out to be always zero. The reason for this is twofold: 


I. First of all, the Uhlmann connection belongs to the 5u(n) Lie algebra, so its trace vanishes and so does the trace 
of its curvature. Hence, using Eq. (61), Chu is zero independently of the associated state pk- This has more to 
do with the mathematical structure of the Uhlmann construction, than with physical properties of pk. 


2. On the other hand, the Chern number is 0 as long as there is a smooth gauge defined along the whole BZ 
[Haiii]. Note that this is the case for the Uhlmann ^(n) gauge, as we can always take the gauge Wk = 
which is well defined provided that pk is not singular at some crystalline momentum k. This is a rather natural 
condition and the usual case in physical situations. 


Additionally, the second condition, also implies that higher order Chern numbers vanish for the Uhlmann connection. 
Thus, Chu = 0 in any case. This result was first proven in Ref. m 

In conclusion, this makes not obvious the extension of 2D topological invariants by means of the Uhlmann approach. 
We will show specifically the way to avoid these obstacles and define a proper topological invariant for general density 
matrices associated to these topological systems. 
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FIG. 5: This 2D toms represents the Brilloiun Zone (BZ) where the crystalline momentum of the system leaves. The figure 
represents the partition of the toms in infinitesimal slices of width dkx, computing the Berry (or Uhlmann) phase associated 
to each slice and integrating of the whole homological cycle in kx, one obtains the Chern number (Eq. (62)) or the Uhlmann 
number (Eq. (|67|) respectively. 


B. Topological Uhlmann numbers 


In this section, we explain how to construct single topological invariants associated to density matrices of 2D 
topological insulators and superconductors. The key tool for that is going to be the topological Uhlmann phase 
previously derived. The 2D construction as explained here was originally proposed in Ref. [37] and also equivalently 
in Ref. [36|. A more detailed derivation and discussion not present in previous works is carried out in what follows. 

In the previous subsection, we showed that the Chern number associated to the Uhlmann connection was always 
zero. However, the fact that Chp becomes identically zero and the associated fiber bundle becomes trivial, does not 
imply that all topological properties of density matrices are trivial. If this assertion were true, we could not claim 
that systems at T = 0 can be topologically ordered, as they are just a particular case of generally mixed density 
matrices. And in fact, we see that systems at T = 0 display topological order in many fashions (non-zero Chern 
numbers, topological edge states, Majorana fermions, magneto-electric effect, ...). The key point in the discussion is 
that Chu is not the only topological invariant that we can construct on a torus out of a density matrix. 

With that in mind, let us analyst the problem at T = 0 from a different perspective. Actually, in the Berry case, 
the Chern number (59) can be rewritten as namEni 


ch = T ^ dk, 


d^Bjkx 

dkr 


(62) 


where = § dkyA^{kx^ky) is the Berry phase along the /c^-nontrivial homological circle of the torus at the 

point kx^ and § dkx denotes the integration along the nontrivial homological circle. See Fig. 

To prove the equality (|6^ from the general definition of the Chern number, one divides the surface integral in Eq. 


(59) into small slices along the kx direction. 


Abz 


d‘^kFxy{k), 


(63) 


where Sbz is an infinitesimal element of area that is limited by one ^-nontrivial homological circle of the torus at the 
point kx and another one at the point kx + dkx^ as shown in Fig. The U(l)-gauge, which may be ill-defined over 
the whole BZ, is always well-defined in a sufficiently small slice, so that the Stokes’ theorem is applicable. Thus, 


27 ^ Jsbz 


- / (fkF,y{k) = 


OSbz 


A^{k)dk = + ^B{kx + dk^) = 


d^B{k, 

dk:r. 


-dkx 


(64) 


Then, in the limit of slices with infinitesimal width, summing the slices at both sides of Eq 
performing the integral and one easily gets Eq. ( [6^ . 

Note that, Eq. (62) gives indeed a lot of information. If ^B{kx) is a continuous function and 4 >b(A:; 
^B{ka 


(64) is equivalent to 


0) = 


27r), then we trivially get Ch = 0. If ^B{kx) displays some 27r-discontinuous jumps along the kx circle, then 
we can always write 




(65) 
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(a) (b) 

FIG. 6: A pictorial image of the QWZ model in real-space is shown here. In (a), we represent the on-site <7z part of the 
Hamiltonian interaction, where the arrows represent the spins in z—direction. In (b), we show an example of the effect of the 
spin-dependent hopping in the interaction. In particular, |t)^ —^ |t)^ when hopping along the x—axis, but It)^ —^ \i)^ 

when hopping along the y—axis. 


where is a smooth function equal to § dkyA^{kx^ky) mod. 27r, k^ is the point at which the discontinuity 

happens and 0 denotes the Heaviside step function. This way we can compute the derivative and perform the 
integration along the nontrivial homological circle of the torus. Actually, the Chern number is measuring the 
number of these 27r jumps, which is equivalent to the number of windings of ^B{kx) as the k^ circle is covered. This 
is indeed a topological invariant. In particular a winding number, that associates every state of the system with an 
homotopy class of the Berry phase mapping 

^ (66) 


between the nontrivial homological circle along kx and the complex Berry phases U{1) = S^. 

As we have seen before, Eq. (59) does not admit a direct generalization to density matrices, however, the equivalent 


formula (62) for the Chern number at T = 0 and pure states allows for a nontrivial generalization to density matrices. 


To that aim, we will just replace the Berry phase ^B{kx) by the Uhlmann phase ^ij{kx)^ Eq. (47), in Eq. (62): 


nu 


^ / dh 

2^7 ■ 


d^ujkx) 

dkr 


(67) 


Note that, many of the properties of the Berry case are kept. This number np is also an integer by the same 
arguments given previously for the Chern number. If the ^\j{kx) is a continuous function along the kx circle and 
^\j{kx ^ 0) = ^\j{kx 27r), then we trivially obtain np = 0. Otherwise, np measures again the number of 27r- 
discontinuous jumps along the nontrivial homological circle. It is also trivially gauge invariant, as the Uhlmann 
phase is gauge invariant already. 

This integer number is therefore a topological invariant which classifies the density matrices of a quantum system 
according to the homotopy class of the Uhlmann phase mapping. 




( 68 ) 


between the nontrivial homological circle along kx and the complex Uhlmann phases U{1) = S^. 

Eurthermore, the limit at zero temperature is also well-defined. Eor pure states, the Uhlmann phase coincides with 
the Berry phase ff>p = Thus, if we compute np in a thermal (Gibbs) state, then the following relation is trivially 
satisfied 


np 


T^O 


> Ch, 


(69) 


hence, the generalization is faithful. 

Additionally, since ff>p is an observable and we even discussed certain general schemes to measure it in different 
platforms [35], then np is also an observable. We will refer to this topological invariant np as the Uhlmann number. 

As was shown in Ref. the Uhlmann number characterizes representative models for topological insulators 

and superconductors in 2D. Remarkably, in models with high Chern number, one finds novel thermal-topological 
transitions between two nontrivial phases solely driven by temperature. 

In what follows, we will analyze a new model for a two-band topological insulators at finite temperature according 
to the topological invariant np. 
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FIG. 7: Topological phase diagram given by the Uhlmann number for the QWZ model. Varying both the energy imbalance m 
and the temperature T we can go from a topological region to a trivial one. Note that nu remains quantized everywhere. If 
we fix m and vary the temperature T, there is always a critical point where nu goes discontinuously and abruptly to zero. 


C. QWZ Chern insulator. 


As an example of the applicability of this method, we present a model for a 2D topological insulator proposed 
in [53l [54]. Although the model is mathematically simple, it contains the main futures of the physics of a Chern 
insulator. 

The model consist of a simple square lattice of spin-1/2 fermions, with spin-dependent nearest-neighbour hoppings, 
and an energy imbalance between spin up and down (see Fig. |^. The Hamiltonian in real space is given by 


H = 






^f-\-X T d 


CTz - ICT, 


-{-h.c. -f m 


ctcTz 


(70) 


where ct, Cf are fermionic creation and annihilation operators and x and y are unitary vectors along the x and 
y direction respectively. The Pauli matrices (7x,y,z act on the spin-1/2 degree of freedom. The model describes 
the quantum anomalous Hall effect. The (Tx and (jy terms appear due to strong spin-orbit coupling and the cTz 
part represents some ferromagnetic polarization. It has been shown to be realisable in Hg]^_^Mna,Te/Cdi_a,Mna,Te 
quantum wells with a proper amount of Mn spin polarization [92| . 

In Fig. [^we sketch the dynamics of the system for a certain spin orientation as an example. Figure |^a) represents 
the on-site cFz part, and Fig. |^b) shows the hopping mechanism for a certain spin state in order to visualise the 
dynamics. When a fermion with spin up in the z—direction hops along the x—direction, then |t)y- If if 

hops along the ^-direction, then |t)^ —as it is shown in Fig. [^b). 

Taking periodic boundary conditions the Hamiltonian in momentum space reads 


H(k) = sin {kx)(Tx sin {ky)ay -h [m-\-cos{kx) cos {ky)]az, 


where m is some energy imbalance between to atoms of different species. 


At T = 0, the different topological phases as classified by the Chern number, Eq. (59), are: 


{ +1, if — 2 < m < 0, 
— 1, if 0 < m < 2, 

0, if otherwise. 


(71) 


(72) 


For nontrivial regions Ch = ±1, the system presents metallic edge states at the boundary. Indeed, the conductivity 

2 

of these edge states is given by the Chern number (lEdge = 

At finite temperature, the different topological phases as classified by the Uhlmann number, Eq. (67), are graphically 
represented Eig. The system displays nontrivial topological phases np = ±1 even at nonzero temperature provided 
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^u(^cc) 



FIG. 8: We represent $u(fc) in the QWZ model for m = —1. On the left figure, we set T = 0.1 and observe a discontinuous 
jump in the phase which is responsible for the robustness of the topological phase at finite temperature. On the right figure, we 
take a higher temperature T = 1.2 where the system is not topological anymore, and the discontinuity in the phase disappears. 


it is below a certain critical value Tc, where np goes to zero. This critical temperature Tc reaches its maximum at 
the points m = ±1 within the topological phases Ch = ±1 at T = 0. These points are associated with the highest 
value of the gap. This matches with the common intuition that the resilience of topological order has to do with the 
ratio between the gap and the temperature, and the larger the gap, the more stable it is. As expected, in the limit of 
T = 0 we recover the same topological diagram as given by the Chern number. 

Actually, we can understand even better the discontinuity in the Uhlmann number when crossing different topo¬ 
logical regions looking at Fig. We represent ^u(^cc) for different topological sectors. Namely, we have set m = —1 
which is a topological region at T = 0. At low temperature we see that there is still a jump in the phase at a certain 
critical momentum kc = 0, which accounts for np = 1. However, in the same figure we can see what happens at higher 
temperature where the jump disappears, the Uhlmann phase is a continuous function of kx and hence, up = 0. This 
supports the previous explanation that the topological Uhlmann number accounts for the number of discontinuous 
jumps of the Uhlmann phase in momentum space. 

Thus, we see that thermal topological phase transitions are not a unique phenomenon of the ID case [35] and they 
are also found in two-dimensional systems. 


V. CONCLUSIONS 

In this paper we have characterized a new emblematic model for a topological insulator, the Qi-Wu-Zhang model 
introduced in [SI EH using the topological Uhlmann number introduced in m in order to find a stable SPT phase 
under external thermal fluctuations in two-dimensions. This was first originally discovered in emblematic examples of 
topological insulators and superconductors in one-dimension m, which later lead to higher dimensional constructions 
ESIET]. A complete phase diagram of the topological Uhlmann number as a function of the parameters of the 
Hamiltonian in the Qi-Wu-Zhang and as a function of the temperature of the system has been shown in Fig. We 
clearly see stable thermal topological phases in large regions of the phase diagram where the topological invariant 
remains quantized and critical temperatures where the system jumps onto a trivial phase with zero value for the 
topological Uhlmann number. 

In order to apply our newly developed topological invariant |37j to the topological insulator in 2D given by the 
Qi-Wu-Zhang model, we have presented a self-contained explanation of the theory of geometric phases in quantum 
mixed states represented by density matrices and how they may account for new topological effects under certain 
conditions. Thus, we have introduced the main concepts and results of the Uhlmann geometric theory. We have used 
the Uhlmann geometric phase as a tool to characterise SPT systems like topological insulators and superconductors at 
finite temperature. The Uhlmann phase acquires a topological character in this situations remaining quantized even 
for a finite range of temperature. In fact, when computed for a thermal state of generic two-band models in ID, the 
complicated expression of the Uhlmann phase simplifies significantly and can be related to the winding number of the 
Hamiltonian and the temperature of the system. It has also a well-defined limit at T = 0, tending to the Berry/Zak 
phase, which makes the theory consistent and faithful. 

Morover, we have presented how this tool can be generalised to 2D systems in a detailed fashion. The topological 
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invariant in this case is the topological Uhlmann number, which generalises the Chern number for symmetry-protected 
time-reversal-broken topological insulators and superconductors in 2D. In fact, with these results, we may think of 
the possibility of extending the current classification of topological insulators and superconductors on several spatial 
dimensions [86l [87] (also called the “periodic table”), to the situation of thermal SPT states with fermions using 
topological Uhlmann numbers. 

Let us emphasise that the existence of critical temperatures seems rather natural in the Uhlmann approach. We 
find that for thermal states, this approach sets on equal footing the temperature and the coupling constants in the 
original Hamiltonian. Therefore, if there is a critical behavior as a function of tunnelings and/or staggered potentials, 
then one should obtain a critical behavior with temperature as well. Moreover, since by raising the temperature T, 
the quantum coherence properties of any state are affected, it is rather natural that the Uhlmann number decreases 
by warming the system up. 

Beyond the scope of this paper it is the possible adaptation of the Uhlmann formalism to explore interacting 
topological phases at finite temperature. They include not only bosonic/spin topological phases but also interacting 
fermionic SPT phases that have been recently proposed and classified [93H95] based on group cohomology and su¬ 
percohomology theory, and response field theory corresponding to a non-trivial fiber bundle of the symmetry group 
[961197]. In this regard, the application of topological Uhlmann invariants to more complicated models and higher 
dimensions may lead to novel thermal topological transitions and new effects. 
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